Adhesion of membranes via actively switched receptors 
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We consider a theoretical model for membranes with adhesive receptors, or stickers, that are 
actively switched between two conformational states. In their 'on'-state, the stickers bind to ligands 
in an apposing membrane, whereas they do not interact with the ligands in their 'off'-state. We 
show that the adhesiveness of the membranes depends sensitively on the rates of the conformational 
switching process. This dependence is reflected in a resonance at intermediate switching rates, which 
can lead to large membrane separations and unbinding. Our results may provide insights into novel 
mechanisms for the controlled adhesion of biological or biomimetic membranes. 

PACS numbers: 87.16.Dg, 05.40.-a, 87.68. +z 



Introduction - Fluid lipid bilayers are the fundamental 
building blocks of biological membranes and, as such, 
have been studied extensively 0. As a consequence, 
many equilibrium properties of lipid membranes are now 
well understood. However, biological membranes are of- 
ten out of equilibrium, driven by embedded active com- 
ponents or by the coupling to the dynamic cytoskeleton. 
These active processes affect the membrane behavior, 
e.g. by enhancing membrane shape fluctuations 0, 0, H 
or lateral diffusion ||. The enhancement of shape fluc- 
tuations has been theoretically predicted for membranes 
with ion pumps as active force centers , and experimen- 
tally verified for giant lipid vesicles with embedded bacte- 
riorhodopsin pumps 3] . Other active membrane proteins 
have been suggested to couple to the local thickness Q 
or curvature p| of the membranes. 

In this letter, we consider another type of active mem- 
brane process. The active components embedded in 
these membranes are adhesion molecules, or 'stickers', 
that can be switched between two states: an 'on'-state 
in which the stickers interact with a second membrane, 
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FIG. 1: A membrane with active stickers (top) adhering to 
an apposing membrane. In one of their two conformations 
(stretched), the stickers are 'on' and can bind to ligands in 
the second membrane. The stickers are actively switched be- 
tween their conformations, with on- and off-rates that are 
independent of their binding state. Switching off a bound 
sticker (left) requires an energy input equivalent to the sticker 
binding energy, which keeps the system out of equilibrium. 



and an 'off'-state in which this interaction is negligible 
(see Fig. ^| . An important biological example of actively 
switched adhesion molecules are integrins. In one of their 
conformations, the integrin molecules are extended and 
can bind to ligands present in an apposing membrane 
In another conformation, the molecules are bent 
and, thus, deactivated. The transitions between these 
conformations are triggered by signaling cascades in the 
cells, which typically require energy input, e.g. via ATP. 
In biomimetic applications, conformational transitions of 
appropriately designed stickers may also be triggered by 
light 0. 

In an equilibrium situation, the adhesiveness of a 
biomembrane only depends on the concentration and 
binding energy of the sticker molecules [lo|. We show 
here that the adhesiveness of membranes with active 
stickers also depends strongly on the switching rates of 
the stickers. 

Model - Before considering extended membranes with 
active stickers, let us first focus on a single sticker. In the 
simple example shown in Fig. |2 the sticker can only at- 
tain two different separations from its ligand. The sticker 
is within binding range of the ligand at close separation, 
and out of binding range at a larger separation. At both 
separations, the sticker can be in either of its two confor- 
mations, 'on' or 'off'. In this one-sticker model, we as- 
sume that the change in elastic energy of the membrane, 
which is associated with the change of separation, is neg- 
ligible. The energy of an 'off '-sticker is then independent 
of the separation, which implies that the rate for the tran- 
sition from the larger to closer separation is identical to 
the transition rate for the reverse process. In the dimen- 
sionless units used in this example, the two transition 
rates are equal to 1. An 'on'-sticker, in contrast, gains 
the dimensionless binding energy u when bound to the 
ligand. If we adopt standard Metropolis dynamics, the 
dimensionless unbinding rate of this sticker is exp(— u), 
whereas the rate for the reverse process is 1. 

Let us now consider the transition rates between the 
'on'- and 'off'-state. The rate uj+ for switching the 
sticker on is independent of the membrane separation 
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if we assume that this transition is thermally activated 
and that the corresponding energy barrier is not affected 
by the presence of the ligand. However, the off-rates 



for a bound sticker and 



for an unbound sticker 



in general can be different. In equilibrium, the transi- 
tion rates between the states have to obey detailed bal- 
ance, which implies that the total rate u> b _uj + for passing 
through the four states in a clockwise direction has to 
be equal to the total rate exp(— u)u/" tv + for the counter- 
clockwise cycle. This leads to the equilibrium condition 
uo b _ = exp(— m)cj" . 

In the following, we will consider a situation in which 
the stickers are switched off by an active process that 
violates detailed balance and, thus, the equilibrium con- 
dition above. The active process requires energy input, 
e.g. via ATP hydrolysis. The energy that is set free by the 
hydrolysis of an ATP molecule typically is significantly 
larger than the binding energy of an adhesion molecule. 
Therefore, it seems not unreasonable to assume that the 
rate of an ATP-driven off-switching process is indepen- 
dent of whether the sticker is bound or not. We then 
have u} b _ = = (jj_.^lj At small ATP concentration, 
the 'off '-rate <x>_ should be proportional to this concen- 
tration. 

Let us now consider a simple, discrete model for a 
membrane with active stickers that interact with ligands 
in an apposing membrane. In this model, a 'reference 
plane' parallel to the membranes is discretized into a 
square lattice with lattice spacing a. At each site i, a 
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FIG. 2: A simple four-state model of an active sticker. In this 
model, the sticker can attain only two different separations 
from its ligand, and two conformational states, 'on' or 'off'. 
At close separation, an 'on '-sticker is bound to its ligand with 
binding energy U — uksT. The transition rates between the 
four states are given in units of 1/r, where r is a typical 
timescale for the local membrane relaxation between the two 
separations. The dimensionless rates uj h _ and uft for switching 
a sticker off may in general depend on whether the sticker is 
bound or unbound, whereas the 'on'-rate u>+ is taken to be 
independent of the separation. 



sticker molecule is present that can be in the state 'on' 
{rii = 1) or 'off' (rii = 0). The sticker distribution within 
the membrane thus is approximated by a square-lattice 
distribution. The conformations of the membranes are 
described by the local membrane separation U > at the 
sticker sites i. The membrane separation is nonnegative 
since the membranes cannot penetrate each other. 

The configurational energy of the membranes at time 
t has the form 

H{l,n(t)} = E [^2 ( A ^) 2 + Mt)V(h)\ (1) 

i 

The first term in this configurational energy represents 
the bending energy of the membranes. Here, n — 
«iK2/(ki + K2) is the effective bending rigidity of the 
membranes with rigidities K\ and K2, and the discretized 
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l x ,y+i + lx,y-i — ^lx,y captures the local mean curvature 
of the separation field I. The second term is the time- 
dependent interaction energy of the membranes. If the 
sticker at site i is 'on' at time t, then rii(t) is 1 and 
the membranes interact with the short-ranged sticker 
potential V(U). For simplicity, we model the short- 
ranged sticker interaction by the square-well potential 
V(h) = -U for < h < l r and V(h) = for h > l r 
with characteristic binding energy U > and potential 
range l r . To reduce the number of parameters, we use 
the rescaled separation field Zi ~ {U/a)^J ' k/ (ksT) in our 
simulations. In the following, the rescaled potential range 
z r = {l r / a)\J k/ '(ksT) of the stickers has the fixed value 
z r = 0.1. 

The dynamics of the model here is simulated via a 
master equation with discrete time steps r. In a time 
step r, the rescaled separation Zi at each lattice site i is 
shifted to a value z\ with standard Metropolis transition 
probabilities: the 'new' value z[ is adopted with proba- 
bility 1 if the 'new' conformation is lower in energy, and 
adopted with probability exp[— (H 1 — i7)/(fcsT)] if the 
energy H 1 of the 'new' conformation is higher than the 
energy H of the 'old' conformation. Below, we choose 
z[ = Zi + C,p where p is a random number between -1 and 
1, and ( = 0.5 is the step size. 

To ensure consistency with versions of the model that 
are continuous in time [l2j , the sticker interaction energy 
is averaged over each time step r. If the sticker was al- 
ways 'on' in the interval r, the interaction energy is V(zi). 
If the sticker was always 'off', the interaction energy is 
0. If it was 'on' during a fraction q of the interval r, the 
interaction energy is qV(zi). The membrane separation 
field thus evolves in a stochastic sticker potential that is 
averaged over the time steps t. 

Results - The attractive membrane interaction medi- 
ated by the 'on'-stickers acts against a steric repulsion 
of the membranes caused by shape fluctuations. If the 
sticker-mediated attraction is weak, the membrane fluc- 
tuations will lead to large membrane separations, and 
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FIG. 3: The averaged rescaled separation z of the membranes 
as a function of the mean switching rate ui for the average frac- 
tion X = 0.5 of 'on'-stickers. The membrane separation has a 
maximum at intermediate switching rates oj, and descreases 
with increasing sticker binding energy U . The data points 
are obtained from Monte Carlo simulations with up to 10 7 
steps per lattice site on lattices with up to 140 x 140 sites. 
Statistical errors are smaller than the symbol sizes. 



eventually to an unbinding of the membranes. The av- 
erage membrane separation in the steady state and the 
location of the unbinding transition therefore reflect the 
adhesiveness of the membranes. Clearly, the adhesiveness 
should increase with the sticker binding energy U and 
the average fraction of 'on'-stickers X = oj+/(lu + + cj_). 
However, even for a constant binding energy and a con- 
stant fraction of 'on'-stickers, the average membrane sep- 
aration depends sensitively on the switching rates. In 
Fig. [3J the average membrane separation in the steady 
state is shown as a function of the mean switching rate 
oj = t>(oj + + ojJ). The fraction of 'on'-sticker is X = 0.5, 
which implies oj + = oj_ . The average membrane separa- 
tion exhibits a resonance: it is maximal at intermediate 
switching rates oj, and considerably smaller at high or low 
switching rates. This resonance of the membrane sepa- 
ration can be understood from an associated resonance 
in the escape rate of a single bound sticker molecule, see 
Appendix. 

The resonance of the membrane separation can lead to 
a reentrant unbinding transition, see Fig. For specific 
values of the fraction X and binding energy U of 'on'- 
stickers, the membranes then first unbind with increasing 
switching rate oj, and later rebind. The unbinding tran- 
sitions of the active membranes are continuous. At an 
unbinding point, the average separation z continuously 
tends to infinity [lflj . and the contact probability P&, the 
average fraction of lattice sites with z\ < z r , continuously 
tends to zero. 

The steady-state phase behavior of the membranes can 
be characterized by unbinding lines, see Fig. El The reso- 



nance at intermediate switching rates is reflected in min- 
ima of the inverse critical binding energy l/U c . This 
resonance effect increases with decreasing fraction X of 
'on'-stickers. In the asymptotic limits of small and large 
switching rates, the unbinding behavior of the active 
membranes can be related to the behavior of 'equilib- 
rium membranes'. In the limit of small switching rates, 
the active membranes are equivalent to membranes with 
a concentration X of stickers with 'permanent' interac- 
tion potential V(z). The effective adhesion potential 
V c { of such membranes can be determined exactly by 
integrating over the stickers degrees of freedom in the 
partition function |lf)l |. Close to the unbinding point, 
this effective potential is given by exp [—V f{zi) / ksT] = 
Xexp [— V{zi)/ksT] + 1 — X. Hence, the critical sticker 
binding energy is 



U c {oj = 0) = -k B T\n- 



X 



exp [Ul im /k B T] + X - 1 



(2) 



where U^ m is the critical binding energy of a homoge- 
neous membrane in a square- well potential with the same 
rescaled potential range z r — 0.1. From Monte Carlo sim- 
ulations, we obtain U^ m = 1.26 ± 0.01. In the limit of 
fast sticker switching, the critical binding energy behaves 
as 



U c {oj = oo) = U^ m /X 



(3) 



since the membranes then move in the average potential 
XV(zi). 

Discussion and conclusions - We have shown that ac- 
tive switching processes of stickers can have a strong im- 
pact on the separation and unbinding behavior of mem- 
branes. The adhesiveness of an 'equilibrium sticker' is 
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FIG. 4: Average membrane separation z (filled circles) and 
contact probability Pb (open circles) as a function of the mean 
switching rate oj. The fraction of 'on'-stickers here is X = 
0.2, and the sticker binding energy is U = lOfcsT. With 
increasing switching rate, the membranes first unbind and 
then rebind at the critical rates logcj c r = —0.79 ± 0.07 and 
logw c T = —0.02 ± 0.05 (dashed lines), which are obtained 
from extrapolation to Pt — and l/z = 0. 
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FIG. 5: Unbinding lines obtained from Monte Carlo simu- 
lations for several values of the fraction X of 'on'-stickers. 
Here, u c = U c /k,BT is the dimensionless critical binding en- 
ergy of the stickers, and u) is the mean switching rate. At a 
given value of X, the membranes are bound for sticker bind- 
ing energies U larger than U c and, thus, for 1/U < 1/U C , 
and unbound for 1/U > 1/U C - The resonance at intermedi- 
ate switching rates is is reflected in small values of 1/U C and 
most pronounced at X = 0.2. The dashed lines indicate the 
asymptotic limits given by the eqs. © and 



characterized by its binding energy U. In contrast, the 
adhesiveness of an active sticker depends also sensitively 
on the transition rates of the conformational switching 
process. 

At low switching rates, a bound sticker can unbind ei- 
ther by (i) 'hopping' out of the potential well, or by (ii) 
diffusing out of the potential range after being switched 
'off', since the probability for being switched 'on' in the 
subsequent time steps is small. The typical time scales 
for these processes are (i) the intrinsic unbinding time 
tub ~ rex p[C^], and (ii) the inverse 'off '-rate 1/ui-. For 
w_ <C l/t u b> the adhesion behavior is nearly unaffected 
by the sticker switching. The impact of the active switch- 
ing becomes clearly noticable for switching rates with 
w_ ~ 1/iub- For biological adhesion molecules such as 
integrins or selectins, the typical dissociation times i u b 
are of the order of seconds |b|. Switching processes on 
second or subsecond timescales therefore should allow to 
control the adhesion behavior of biological or biomimetic 
membranes. 

For a given concentration X of 'on'-stickers, the 
adhesion strength in our model is minimal at resonant 
mean switching rates u ~ 1/r, see Figs. 01 and 03 The 
characteristic timescale r corresponds to the typical 
relaxation time of membrane segments with linear size 
a, the distance between the stickers. If hydrodynamic 
damping dominates, this relaxation time scales as 
where rj is the viscosity of the surround- 
ing solvent. For the typical bending rigidity k = lOfc^T, 
an average lateral sticker distance of 180 nm, and the 



typical value r/ = 0.06 poise for the cytoplam viscosity 
of a cell, we obtain a relaxation time scale of 0.1 ms. 
On this length scale, intermonolayer friction leads to the 
same estimate for the relaxation time [lfij ]. 

B.R. thanks Marek Napiorkowski for stimulating dis- 
cussions. 
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APPENDIX 

In our simulations, the interaction energy for a Monte 
Carlo step at time t and site i with rescaled separation 
Zi is (V(zi)/T) f t n{t')dt'. This interaction energy in- 
cludes an average over the sticker switching process n(t ) 
in the time interval from t — r to t. To understand the 
resonance in the membrane separation observed in our 
simulations, it is instructive to consider the escape of a 
single sticker, or single particle, from the dimensionless 
stochastic potential v — (u/t) L n(t')dt' with binding 
energy u = U /ksT and < v < oo. At t = 0, the sticker 
is 'on'. To determine the average escape rate k = (e~ v ) at 
time t = t, we consider the stochastic differential equa- 
tion dv/dt = (zt/r)n(t) and construct the corresponding 
Fokker-Planck equations [17] 



5P ° P P 



p p u dPl 



(4) 
(5) 



Here, P\{v,t) and Po(v,t) are the probability distribu- 
tions for n(t) — 1 and n(t) — 0, respectively. Multiplying 
each of these equations with e~ v and integrating over v 
leads to the ordinary first-order differential equations 

bi = uj + bo — (o;_ + u/r)bi , bo — — uj + b$ (6) 

for the two auxiliary variables b\(t) — J e~ v Pi(v, t)dv 
and bo(t) — J e~ v Po(v, t)dv. In terms of these auxil- 
iary variables, the average escape rate is k = (e~ v ) = 
&i(t) + &o(t). For the initial condition n(0) = 1, and 
thus 6i(0) = 1 and bo(0) = 0, we obtain 



cosh A 



2ojt 



2A 





U + 2(jJT 


j exp 


2 



(7) 



with A 



For the fraction X - 
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0.5, the escape rate k as a 
function of lj is shown in Fig. 6. In the asymptotic 



limit of slow switching or small u>, the sticker has a 
high probability to be still in the initial 'on'-state. 
The average potential of the sticker then is identical 
with its binding energy u, and the escape rate simply 
is e~ u . In the limit of fast switching or large u>, the 
sticker will switch very often between the states. The 
average potential then is Xu = u/2, and the escape rate 
e~ Xu ~ e~ w l 2 . But at the intermediate switching rates, 
the average potential depends on the switching process 
and can vary from values close to zero (if the sticker is 
switched 'off' early in the interval r and remains 'off') 
to the value u (if the sticker remains in the 'on'-state). 
The resonance effect at the intermediate switching rates 
can be understood from the contribution of switching 
processes with small average potential v to the escape 
rate k — (e~ v ). The effect is closely related to the 
'resonant activation' of particles in potentials with a 
fluctuating barrier [18]. 
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FIG. 6: Escape rate k for a single sticker as a function of the 
mean switching rate u> at the average fraction X — 0.5 of the 
'on'-state and several values of the sticker binding energy U , 
as given by eq. @- 



